Decoherence of Quantum- Enhanced Timing Accuracy 
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Quantum enhancement of optical pulse timing accuracy is investigated in the Heisenberg picture. Effects of 
optical loss, group-velocity dispersion, and Kerr nonlinearity on the position and momentum of an optical pulse 
are studied via Heisenberg equations of motion. Using the developed formalism, the impact of decoherence by 
optical loss on the use of adiabatic soliton control for beating the timing standard quantum limit [Tsang, Phys. 
Rev. Lett. 97, 023902 (2006)] is analyzed theoretically and numerically. The analysis shows that an appreciable 
enhancement can be achieved using current technology, despite an increase in timing jitter mainly due to the 
Gordon-Haus effect. The decoherence effect of optical loss on the transmission of quantum-enhanced timing 
information is also studied, in order to identify situations in which the enhancement is able to survive. 

PACS numbers: 42.50.Dv, 42.65.Tg, 42.81.Dp 



I. INTRODUCTION 

It has been suggested that the use of correlated photons is 
able to enhance the position accuracy of an optical pulse be- 
yond the standard quantum limit, and the enhancement can be 
useful for positioning and clock synchronization applications 
Ql]. Generation of two photons with the requisite correlation 
has been demonstrated experimentally by Kuzucu et al. [2], 
but in practice it is more desirable to produce as many corre- 
lated photons as possible in order to obtain a higher accuracy. 
To achieve quantum enhancement for a large number of pho- 
tons, a scheme of adiabatically manipulating optical fiber soli- 
tons has recently been proposed [3], opening up a viable route 
of applying quantum enhancement to practical situations. The 
analysis in Ref. [3] assumes that the optical fibers are loss- 
less, so the Heisenberg limit [4] can be reached in principle. 
In reality, however, the quantum noise associated with optical 
loss increases the soliton timing jitter and limits the achievable 
enhancement. Compared with the use of solitons for quadra- 
ture squeezing [5], the adiabatic soliton control scheme poten- 
tially suffers more from decoherence, because the soliton must 
propagate for a longer distance to satisfy the adiabatic approx- 
imation. The effect of loss on a similar scheme of soliton mo- 
mentum squeezing has been studied by Fini and Hagelstein 
01, although they did not study the timing jitter evolution rel- 
evant to the scheme in Ref. [3], and did not take into account 
possible departure from the adiabatic approximation. 

In this paper, the decoherence effect of optical loss on the 
timing accuracy enhancement scheme proposed in Ref. ||3J is 
investigated in depth, in order to evaluate the performance of 
the scheme in practice. Instead of approachingthe problem in 
the Schrodinger picture like prior work 01 H, 0,01, tn ' s P a P er 
primarily utilizes Heisenberg equations of motion, since they 
are able to account for dissipation and fluctuation in a more el- 
egant way. For simplicity, scalar solitons, as opposed to vector 
solitons studied in Ref. yfl, are considered here. The theoret- 
ical and numerical analyses show that, despite an increase in 
timing jitter due to quantum noise and deviation from the adi- 
abatic approximation, an appreciable enhancement can still be 
achieved using a realistic setup. 

The developed formalism is also used to study the propa- 
gation of an optical pulse with quantum-enhanced timing ac- 



curacy in a lossy, dispersive, and nonlinear medium, such as 
an optical fiber, in order to identify situations in which the 
enhancement can still survive. The effect of loss on many 
correlated photons sent in as many channels has been inves- 
tigated by Giovannetti et al. (U], but their analysis focuses on 
a relatively small number of correlated photons and does not 
include the effects of dispersion and nonlinearity. 

This paper is organized as follows: Section HI1 defines the 
general theoretical framework, and derives the standard quan- 
tum limits and Heisenberg limits on the variances of the pulse 
position and momentum operators. Section[nI]studies the evo- 
lution of such operators in the presence of loss, group-velocity 
dispersion, and Kerr nonlinearity, and determines the effect of 
dissipation and fluctuation on the position and momentum un- 
certainties. Section [IV] investigates theoretically and numeri- 
cally the impact of optical loss on the adiabatic soliton con- 
trol scheme using realistic parameters, while Sec. [V] studies 
the decoherence effect on the transmission of the quantum- 
enhanced timing information in various linear and nonlinear 
systems. 



II. THEORETICAL FRAMEWORK 

A. Definition of pulse position and momentum operators 

The positive-frequency electric field of a waveguide mode 
at a certain longitudinal position can be defined as |§] 
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where n is the refractive index, 77 is the real part of n, S is the 
transverse area of the waveguide mode, and c(a>) is the photon 
annihilation operator. The annihilation operator is related to 
the corresponding creation operator via the commutator |8[], 



[c(co),c + (co')] = 5(o)-C(j')- 



(2) 



For a pulse with a slowly-varying envelope compared with the 
optical frequency, the coefficient in front of the annihilation 
operator can be assumed to be independent of frequency and 
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can be evaluated at the carrier frequency C0q, so that the elec- 
tric field is proportional to the temporal envelope annihilation 
operator A (t), 
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1 



A(t) = 
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d(co) = c((0 + a>o) 
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(3) 
(4) 
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The temporal envelope operator A (f) and the spectral operator 
fl(ffl) evidently also satisfy the following commutation rela- 
tions with their corresponding creation operators, 



[A(t),AHt')} = 8(t-t'), 
[d((o),d f (co')] = 8(co - co'). 

The total photon number operator can be defined as 



N = J dtA f (t)A(t), 
and the pulse center position operator as ^ 



where 



T = - I dttA r (t)A(t), 



N=(N) 
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is the average photon number. This definition uses 1 /N as 
the normalization coefficient, instead of the inverse photon 
number operator N~ l used by Lai and Haus IflOfl . in order 
to express the position operator in terms of normally ordered 
optical field operators that are easier to handle, as well as to 
avoid the potential problem of applying N~ 1 on the vacuum 
state. As long as the photon-number fluctuation is small, the 
position operator naturally corresponds to the measurement of 
the center position of the pulse intensity profile. An average 
longitudinal momentum operator can be similarly defined, 



Q. = J dco cod\co)d(co) 



(11) 



If the quantum state is close to a large-photon-number coher- 
ent state, A can be approximated as (A) + 8A, with 0(8A) <C 
0(A). Equations (|9]l and (fTTt then become the approximate 
position and momentum operators defined by Haus and Lai 
for soli tons in a linearized approach Illlll . The linearized ex- 
pressions also describe how they can be accurately measured 
in practice using balanced homodyne detection Illlll . 

For simplicity, we shall hereafter assume that (f) = and 
(Q.) = H. In the systems considered in this paper, these 
two quantities remain constant throughout propagation, if t 
is regarded as the retarded time in the moving frame of the 
optical pulse. 



The commutator between the position and momentum op- 
erators is 



[f,Ci] 



iN 



By the Heisenberg uncertainty principle, 
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B. Derivation of standard quantum limits 

The standard quantum limits and Heisenberg limits on 
(A 2 ) and (r 2 ) should be expressed in terms of the pulse 
width Af , defined as 



A; : ( - / dtt 2 A\t).\u 
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(14) 



and the bandwidth Ao, 
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(15) 



To calculate the standard quantum limit on the position uncer- 
tainty, consider the expansion 

(£2 2 ) = ^(^Jd(o (Od^d J da' (o'd'^d^j , (16) 

where we have written d = d(co) and d' — d(co') as shorthands. 
Rearranging the operators, 



(£2 2 ) = -(- [dco CO 1 ah 



N \N 

I l,m / dco' cooy'iYu'ua 



(17) 



The first term on the right-hand side of Eq. ( fTTb is proportional 
to A© 2 , while the second term contains a normally ordered 
cross-spectral density. To derive the standard quantum limit, 
we shall assume that the cross-spectral density satisfies the 
factorization condition: 



(18) 



This condition is always satisfied by any pure or mixed state 
with only one excited optical mode, such as a coherent state 
lfl2[ [l3ll . The second term on the right-hand side of Eq. ( fTTb 
becomes 

^ j dco J dco' coco' (d^ d n ad') °c (n) 2 , (19) 
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which is assumed to be zero, as per the convention of this 
paper. Thus, the variance of Cl is 



Aco 2 
~N~ 



(20) 



where the subscript "coh" denotes statistics of coherent fields 
lfl2[[l3Tl given by Eq. dTsb . By virtue of the Heisenberg uncer- 
tainty principle given by Eq. ( fT3b . the standard quantum limit 
on the position variance is hence 



(f 2 ) > (f 2 ) 



SQL 
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(21) 



This limit is applicable to any pure or mixed state, and is con- 
sistent with the one suggested by Giovannetti et al. for Fock 
states [ 1]. A very similar derivation of the limit for Fock states 
and coherent states is also performed by Vaughan et al. |9J]. 

Owing to Fourier duality of position and momentum in the 
slowly-varying envelope regime, the standard quantum limit 
on the momentum can be derived in the same way. The vari- 
ance of T, assuming coherent-field statistics, is 



At 2 



)coh AT 



(22) 



'coh 

and the standard quantum limit on the momentum variance is 

1 1 



(6 2 >SOL = 
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states, and acceptable for any quantum state with a small 
photon-number fluctuation, such as a large-photon-number 
coherent state. We then obtain the following approximate in- 
equality, 



(Cl 2 ) < Aco 2 



(27) 



With the Heisenberg uncertainty principle given by Eq. ( fT3l ) 
and the upper bound on (Cl 2 ) given by Eq. d27l >. one can then 
obtain the Heisenberg limit on the uncertainty of T: 



(f 2 )>(f 2 ) H = 



1 



AN 2 Aco 2 ' 



(28) 



Equation d28l) is again consistent with the Heisenberg limit 
suggested by Giovannetti et al. Jj]], although the derivation 
here shows that it is not only valid for Fock states but also 
correct to the first order for any quantum state with a small 
photon-number fluctuation. 

The Heisenberg limit on (Cl 2 ) is similar, 



1 



AN 2 At 2 ' 



(29) 



A more exact derivation of the Heisenberg limits is given in 
Appendix lAl where the inverse photon-number operator N~ 1 
is used instead of 1 /N in the definitions of T, Cl, Aco, and 
At. The difference between the approximate Heisenberg limits 
derived here and the exact Heisenberg limits in AppendixlAlis 
negligible for small photon-number fluctuations. 



C. Derivation of Heisenberg limits 

To derive the Heisenberg limit on the position uncertainty, 
one needs an absolute upper bound on the momentum uncer- 
tainty (Cl 2 ). Consider the following non-negative quantity 
proportional to the coherence bandwidth squared, 

Jdcoj dco' (CO - CO 1 ) 2 (d^d'^dd') > 0. (24) 

This quantity is non-negative because (co — CO 1 ) 2 is non- 
negative and (d^d'^dd') is also non-negative [13]. It can be 
rewritten as 

J dco J dco' (CO - CO 1 ) 2 (d^d'^da) 
= J dco J dco' (co — co') 2 (a 1 dd n a') , (25) 
and expanded as 

l^j dco j dco' (co 2 + CO 12 -2coco')d^ dd' ] a! J > 0, 

i(n j dcoc» 2 a^-2N 2 (Cl 2 )>0. (26) 

Here we shall approximate with Af, and neglect any photon- 
number fluctuation. This approximation is exact for Fock 



III. OPTICAL PULSE PROPAGATION IN THE 
HEISENBERG PICTURE 

The classical nonlinear Schrodinger equation that describes 
the propagation of pulses in a lossy, dispersive, and nonlinear 
medium, such as an optical fiber, is given by IU4ll 

dA pd 2 A , l2 „ /a 

•aF = b?-*W aA -T A ' (30) 

where t is the retarded time coordinate in the frame of the 
moving pulse, /3 is the group-velocity dispersion coefficient, 
K is the normalized Kerr coefficient, and a is the loss coeffi- 
cient, all of which may depend on z. The phenomenological 
quantized version that preserves the commutator between A 
andA+isJl 

dA Bd 2 A - t -~ /a. 

i Tz = ^W- KA ' AA ~Y A+is - (31) 

A = A(z,t) is the pulse envelope annihilation operator in the 
Heisenberg picture, and £ is the Langevin noise operator, sat- 
isfying the commutation relation 

[s(z,f),s" + (z' \t')] = <x8(z- z')8(t -t'). (32) 

Rewriting the position and momentum operators in Eqs. (O 
and ( fTTT ) in the Heisenberg picture as T(z) and Cl(z) in terms of 
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A(z,t), differenting them with respect to z, and using Eq. OTb . 
their equations of motion can be derived, 



df(z) 

dz 
dd{z) 

dz 



= p(z)Cl(z)+S T (z) 
= Sa(z), 



(33) 
(34) 



where St and §a are position and momentum noise operators 
defined as 



§t(z) 
Sa(z) 



1 



N(z) 
1 



dt tS^(z,t)A(z,t)+H. c, 



(35) 



N[ _ u dts\z,t)^ij- t )A(z.t He. (36) 



and H. c. denotes Hermitian conjugate. If the noise reservoir 
is assumed to be in the vacuum state, the noise operators have 
the following statistical properties, as shown in Appendix iBl 

(S T (z))=0, (§ n (z))=0, (37) 
a(z)At 2 (z) 



(St(z)S t (z')) 
{S T (z)Sn(z') + Sn(z)S T (z / )) 



N(z) 
a(z)Aco 2 (z) 

N(z) 
«(z)C(z) 



Siz-z 1 ), (38) 
5(z-A (39) 



N(z) 



S(z-z'), (40) 



where C(z) is the pulse chirp factor, defined as 



C(z) 



1 



N(z) 



dt A f (z,f) 



A(z,0 



(41) 



The average position {f(z}} and average momentum (£2(z)) 
are constant and assumed to be zero throughout propagation. 
The variance of Q. is then 



<n 2 (z)) = (n 2 (o))+ dz' 



_,a(zOA©V) 
iV(z') 



(42) 



while the variance of f is more complicated due to the pres- 
ence of dispersion, 

(f 2 (z)) = (f 2 (0)) + (t (0)6(0) + 6(0)? (0)) jV/3(z') 



(6 2 (0)) 



n 2 

rfzW) 



,, a(z')Af 2 (z') 
o JV(z') 

o Jo N(z") 



2 / dz^z') / <fe"j3(z") 
Jo JO 

f a(z"')Aco 2 (z"') 

Jo N(z">) ' 



(43) 



Equation (l43l is the central result of this paper. It is similar 
to that derived by Haus for optical solitons using a linearized 
approach flal . but Eq. (03]) is valid for arbitrary loss, arbi- 
trary dispersion profile j3(z), and arbitrary evolution of pulse 
width Af(z), chirp C(z), and bandwidth A<a(z), so that it is 
able to describe the effect of loss on the quantum enhance- 
ment scheme proposed in Ref. [3]. The first term on the right- 
hand side of Eq. (l43b is the initial quantum fluctuation, while 
the second and third term on the right-hand side describe the 
quantum dispersion effect [16]. In an ideal scenario described 
in Ref. |3], (t 2 (z)) remains constant if the net dispersion 
fedz'fiiz') is zero and quantum dispersion is compensated. 
With loss, however, noise introduces a diffusive jitter given 
by the fourth term on the right-hand side of Eq. 1431 . 



(T 2 (z)) D 



dz 



,a (z')At 2 (z') 
N(z') ' 



(44) 



a less well-known chirp-induced jitter given by the fifth term, 



- [ l ,Ja f ^ r',,a(z")C(z") 



<f 2 (z)) c ^ / dz'P(z') dz 



o N{z") 



(45) 



and also the Gordon-Haus timing jitter lfl7ll given by the sixth 
term, 



{f 2 (z)) GH ^2j Q dz'l3(z')J Q dz"P(z") 

t a{z"')A(o 2 {z"') 
N{z'") 



x / dz 
Jo 



(46) 



In most cases considered here, N 3> 1, (f 2 ) At 2 , 
and (n 2 ) <C A© 2 , so one can use the classical nonlinear 
Schrodinger equation, Eq. 1301 . to predict the evolution of 
At(z), C(z), and Ao(z) accurately. The evolution of (f 2 (z)) 
can subsequently be calculated analytically or numerically us- 
ing Eq. d43l and the classical evolution of At(z), C (z), and 
Aco(z), analogous to the linearized approach IUll[l5ll . 

It is worth noting that the chirp-inducedjitter, Eq. 1451 . de- 
pends on the cross-correlation between the position and mo- 
mentum noise in Eq. 1401 . so it can be positive as well as neg- 
ative, but the sum of the three sources of jitter must obviously 
remain positive. 



IV. EFFECT OF LOSS ON ADIABATIC SOLITON 
CONTROL 

A. Review of the ideal case 

Consider the scheme proposed in Ref. [3] and depicted in 
Fig. Q] Assume that the dispersion coefficient of the first fiber 
/3(z) is negative and its magnitude increases along the fiber 
slowly compared with the soliton period. The classical soliton 
solution of Eq. (l30l l. assuming adiabatic change in parameters 
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z= 



0(z) <o, ia( Z )|t 



1 z 3 ^ 



: = L 



z=L+L' 



FIG. 1 : (Color online) Schematic (not-to-scale) of the adiabatic soli- 
ton control scheme. An optical pulse is coupled into a dispersion- 
increasing fiber of length L with a negative dispersion coefficient /3, 
followed by a much shorter dispersion-compensating fiber of length 
L! with a positive dispersion coefficient j3'. 



/3(z)andJV(z),is01 



A(z,t) = A (z) sech 



exp 



IK 



2 Jo 



dz'\A (z!)\ 2 



2T( z y 



kN(z) 



The adiabatic approximation is satisfied when 



< A, 



N(z) 



dp(z)/dz 
where A is the soliton period, 



dN(z)/dz 



-«A, 
a 



n t 2 (z) 
MZ)= 2\W)\ 



(47) 
(48) 

(49) 

(50) 



The root-mean-square pulse width Af(z) and bandwidth 
Aco(z) then become 



Af(z) 



Ato(z) = 



2a/3 
1 



<z) 



V3- 



y/3 KN(z) ' 
1 kN(z) 
2V3\P(z)\- 



(51) 
(52) 



The bandwidth Ao)(z) is thus reduced in the first fiber. If the 
second fiber has a positive dispersion coefficient /3' so that 
the net dispersion is zero (£ dzfi (z) +P'L'= 0), the quantum 
dispersion effect given by the second and third term on the 
right-hand side of Eq. (143) can be eliminated. Furthermore, if 
/3' has a very large magnitude compared with /3 (z) so that the 
second fiber can be very short compared with the first fiber, the 
effective nonlinearity experienced by the pulse in the second 
fiber can be neglected, and A©(z) remains essentially constant 
in the second fiber. In the lossless case, the final timing jitter 
(T 2 (L + L')} is therefore the same as the input (f 2 (0)), but 
Acd(L + L') has been reduced and the standard quantum limit 
on (T 2 (L + L')}, Eq. d21) . is raised. Provided that the initial 
timing jitter of a laser pulse obeys the coherent-field statistics 
given by Eq. (22) , the final timing jitter is 



(f 2 (L + L')) = (f 2 (0)) = 



Af 2 (0) _ 7t 2 /3 2 (0) 
N ~ T k 2 N 3 



(53) 



while the final standard quantum limit is 

1 3/3 2 (L) 



(f 2 (L + L>)) 



S( 3 L 4NAca 2 (L + L>) k 2 N 3 



(54) 



A timing jitter squeezing ratio, analagous to the squeezing ra- 
tio defined by Haus and Lai Hill , can be defined as 



R = 



(t 2 (L + U)) _7r 2 /3 2 (0) 
(T 2 (L + L<)) SQL ~ 9 P 2 (Ly 



(55) 



The factor of n 2 /9 arises because the initial jitter for a sech 
pulse shape is slightly higher than the standard quantum limit 
given by Eq. ( |2"TT i in terms of the bandwidth. As long as /3 (L) 
at the end of the first fiber is significantly larger than the initial 
value, the timing jitter becomes lower than the raised standard 
quantum limit, R becomes smaller than 1, and quantum en- 
hancement of position accuracy is accomplished. This semi- 
classical analysis is valid in all practical cases, where N ^> I, 
R > l/N, (t 2 ) < Af 2 , (Cl 2 ) < Aco 2 , and is consistent with 
the analysis of exact quantum soliton theory in Ref. J3l. R 
is related to the quantum enhancement factor y defined in 
Ref. JH] by R = 1 /y 2 . The semiclassical analysis is no longer 
valid when R is close to the Heisenberg limit l/N, but as the 
next sections will show, owing to decoherence effects, it is ex- 
tremely difficult for the enhancement to get anywhere close to 
the Heisenberg limit. 



B. Numerical analysis of a realistic case 

To investigate the impact of noise and the validity of the 
adiabatic approximation in practice, a numerical evaluation of 
At (z), C(z), Aco(z), and (f 2 (z)), using Eqs. (H and (03) and 
realistic parameters, is necessary. j3(z) is assumed to have the 
following profile used in Ref. 11911 . 



Hz) 



-12.75 ps 2 /km 

l + (L-z)/V 



(56) 



Lp = 1 km is used here instead of the Lp = 1/12 km used 
in Ref. 11911 . in order to satisfy the adiabatic approximation 
for a longer pulse in this example. Other fiber parameters are 
a = 0.4 dB/km, n 2 = 2.6 x 10" 16 cm 2 AV, A eS = 30 ^m 2 HI , 
Ao = 1550 nm, (Oq — 2tic/Xq, so that k = hcoo(coon2/cA s ff). 
L is assumed to be 2 km. A dispersion-compensating fiber 
with j8' = 127.5 ps 2 /km^= 0.4 dB/km, n 2 = 2.7 x 10~ 16 
cm 2 AV, A e ff = 15 jurn 2 Coll. andZ/ = 110 m is used in the 
numerical analysis as the second fiber. The classical nonlin- 
ear Schrodinger equation, Eq. d30) , is numerically solved us- 
ing the Fourier split-step method II14I1 . An initial sech soliton 
pulse with t(0) = 1 ps, N(0) = 1 .9 x 10 7 , and an initial energy 
of 2.4 pj is assumed. 

Figure|2]plots the numerical evolution of pulse intensity and 
spectrum in the two fibers. As expected, the bandwidth is nar- 
rowed in the first fiber and remains approximately constant 
in the second (z > 2000 m), owing to the latter's relative short 
length. Figure[3]plots the evolution of pulse width Af(z), chirp 
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FIG. 2: (Color online) Numerical evolution of pulse intensity (top) 
and spectrum (bottom). The denser plots for z > 2000 m indicate 
pulse propagation in the second fiber. The color codes are in the 
same arbitrary units as the heights of the plots. 
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C(z), and bandwidth A©(z), compared with the adiabatic ap- 
proximation, Eqs. d5TT > and d52l ). The adiabatic approximation 
is evidently not exact, and the pulse acquires a chirp due to 
excess dispersion in the first fiber, leading to slight refocusing 
in the second fiber. The bandwidth is reduced by a factor of 
2.2, as opposed to the ideal factor of 3.6. 

Figure |4] plots the evolution of the diffusive jitter given by 
Eq. d44i >. the chirp-induced jitter given by Eq. d45l >, and the 
Gordon-Haus jitter given by Eq. d46b . It can be seen that al- 
though the Gordon-Haus jitter increases much more quickly 
than the other jitter components in the first fiber, the for- 
mer drops abruptly in the second fiber (z > 2000 m) due 
to the opposite dispersion. This kind of Gordon-Haus jit- 
ter reduction by dispersion management is well known 12111 . 
The chirp-induced jitter component drops below zero in the 
second fiber, but as noted before, the total noise jitter re- 
mains positive. The final jitter values are numerically deter- 
mined to be (f 2 (L + L')) D = 0.71 (f 2 (0)), (t 2 {L + L')) c = 

-0.93 (f 2 (0)), and (f 2 (L + L')) GH = l-42(f 2 (0)), result- 
ing in a total jitter of 



(f 2 (L + L')) = (f 2 (0)) + (f 2 (L + L')) D 



+ (f 2 (L + L') 
= 2.19(f 2 (0)). 

The final squeezing ratio is hence 



(t 2 (L + L>)) 



GH 



(57) 



FIG. 3: (Color online) Evolution of pulse width At(z) (top), chirp 
C(z) (center), and bandwidth Aco (bottom), compared with the adia- 
batic approximation (dash). Plots of At and Aco are normalized with 
respect to their initial values, respectively. 



R = 



(T 2 (L + L>)) 
(f 2 (L + L>)) SQL 

7t 2 (f 2 (L + L')) N{L + L') Aco 2 (L + L') 
' 9 (f 2 (0)) JV(0) Aw 2 (0) 
0.42 = -3.8 dB. 



(58) 



Despite taking into account the increased timing jitter and the 
non-ideal bandwidth narrowing, a squeezing ratio of —3.8 dB 
is predicted by the numerical analysis, suggesting that one 
should be able to observe the quantum enhancement experi- 
mentally using current technology. 



C. Potential improvements 

As shown in the previous section, the Gordon-Haus effect 
contributes the largest amount of noise in the soliton control 
scheme, despite its partial reduction by dispersion manage- 
ment. Its magnitude at the end of the first fiber can be esti- 
mated roughly as 

2 

(ccL). (59) 



<f 2 (L)) 



GH 



( f2 ^))sQL 

As the length of the first fiber must be at least a few times 
longer than the soliton period A for the adiabatic approxima- 
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FIG. 4: (Color online) Evolution of diffusive jitter (top), chirp- 
induced jitter (center), and Gordon-Haus jitter (bottom). All plots 
are normalized with respect to the initial jitter (r 2 (0)). 



tion to hold and for the bandwidth to be significantly reduced, 
L/A is approximately fixed, and the Gordon-Haus jitter can 
be reduced only if a figure of merit, 



FOM = 



1 

aA 



IB. 

% OCT 2 



(60) 



is enhanced. Since this is a rough order-of-magnitude esti- 
mate, a representative value of A, say at z = L, can be used. 
The figure of merit suggests that the performance of the soli- 
ton control scheme can be improved by reducing the pulse 
width, increasing the overall dispersion coefficient, or reduc- 
ing the loss coefficient. 

Reducing the pulse width is the most convenient way of 
obtaining better enhancement, as the adiabatic bandwidth re- 
duction can be achieved over a shorter distance with less loss 
of photons. For example, using t(0) = 500 fs, L = 1 km, 
Lp = 0.3 km, L' = 44 m, and otherwise the same parame- 
ters as in Sec. IIV Bl the squeezing ratio becomes —6.0 dB, 
while using t(0) = 200 fs, L = 500 m, Lp = 1/12 km, and 
L' = 16.2 m gives a squeezing ratio of —7.3 dB. The shorter 
pulse width, however, significantly enhances higher-order dis- 
persive and nonlinear effects. Raman scattering, in particular, 
contributes additional quantum noise because of coupling to 
optical phonons [22]. It is beyond the scope of this paper to 
investigate these higher-order effects, so a more conservative 
pulse width of 1 ps is used in the preceding section. A larger 



overall dispersion coefficient, on the other hand, means that 
more photons or a higher nonlinearity are required for a soli- 
ton to form, so the Raman effect may also become more sig- 
nificant with a larger dispersion coefficient. The Raman effect 
can be reduced by cooling the fiber and reducing the number 
of thermal phonons [22], if it becomes a significant problem. 

Further advance in optical fiber technology should be able 
to increase the figure of merit by reducing loss, since the spe- 
cialty fibers assumed in Sec. IIV Bl have a higher loss than usual 
transmission fibers by a factor of two. Using a = 0.2 dB/km 
instead of 0.4 dB/km in Sec. IIV B| for instance, reduces the 
squeezing ratio to —4.7 dB. Spectral filtering or frequency- 
dependent gain 1 23] provides another way of controlling the 
Gordon-Haus effect, although it adds another level of com- 
plexity to the experimental setup, and it is beyond the scope 
of this paper to investigate how the frequency-dependent dissi- 
pation or amplification might help the quantum enhancement 
scheme. Finally, the design of the setup assumed in Sec. IIV Bl 
is not fully optimized, and further optimization of parame- 
ters, fiber dispersion profiles, and bandwidth narrowing strat- 
egy should be able to improve the enhancement. 



V. EFFECT OF LOSS ON THE TRANSMISSION OF 
QUANTUM-ENHANCED TIMING INFORMATION 

Provided that quantum enhancement of pulse position ac- 
curacy can be achieved, the information still needs to be 
transmitted through unavoidably lossy channels. It is hence 
an important question to ask how loss affects the quantum- 
enhanced information in optical information transmission sys- 
tems. Equation d43l governs the general evolution of the tim- 
ing jitter under the effects of loss, dispersion, and nonlinearity, 
but in order to estimate the relative magnitude of the decoher- 
ence effects and gain more insight into the decoherence pro- 
cesses, in this section Eq. ( f4~3b is explicitly solved for various 
systems and compared with the standard quantum limit. 



A. Linear non-dispersive systems 

Without dispersion, the timing jitter increases only due to 
the diffusive component (f 2 (z)) D - An analytic expression for 
(f 2 (z)) can then be derived from Eq. d43l . as At(z) andAo(z) 
remain constant, 



(f 2 (z))=(f 2 (0)>- 



Af 2 (0) 

W) 



(l-< 



(61) 



If the initial variance obeys coherent-field statistics, that is, 
(f 2 (0)) = Af 2 (0) /N(0) according to Eq. <|22>, the subsequent 
jitter is 



(t 2 (z)) 



coh 



Af 2 (0) 

W 



(62) 



and obeys the same coherent-field statistics but for the re- 
duced photon number N(z). This is consistent with intu- 
ition. On the other hand, in the high loss limit (az S> 1), the 
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term At 2 (0)/N(z) is likely to dominate over the initial jitter 
(t 2 (0)), so in most cases the position of a significantly atten- 
uated pulse relaxes to coherent-field statistics independent of 
its initial fluctuation. This justifies the assumption in Sec.lIVI 
that a laser pulse exiting a laser cavity has such statistics, re- 
gardless of the quantum properties of the pulse inside the cav- 
ity. 

Equation doTT ) can be renormalized as 



(T 2 (z)) SQL 

: R(0)e- az +4Af 2 (0)Aco 2 (0)(l - e - az ). (63) 



When (f 2 ) < Af 2 and (£2 2 ) < Aco 2 , classical theory pre- 
dicts that 4Af 2 A(B 2 ps 1. Equation d63l then suggests that the 
relative increase in timing jitter is independent of the initial 
squeezing ratio R(0). This is nevertheless not true in general, 
as Af may depend on both Aco and R when the classical theory 
fails. In AppendixICl the exact dependence of Af on Aco and R 
is calculated for a specific multiphoton state with a Gaussian 
pulse shape called the jointly Guassian state. The expression 
4Af 2 Aco 2 is given by 



a a 2* 2 R (1- W 

4ArAor = — 1 



iV 1-1/(JV/?)' 



(64) 



which results in the following exact expression for an initial 
jointly Gaussian state, 



R(z) = R{0)e- az 



R (1-1/AO 2 



N l-l/(NR) 



(l-e- az ). (65) 



For a large photon number (N 3> 1) and moderate enhance- 
ment (1 > R > 1 /AO, 4Af 2 Ao 2 w 1, as classical theory would 
predict for a Gaussian pulse. In this regime, the quantum- 
enhanced information is just as sensitive to loss as standard- 
quantum-limited information. When R gets close to the 
Heisenberg limit 1 /N, however, Af A© approaches infinity. 
This is because maximal coincident-frequency correlations 
are required to achieve the Heisenberg limit [1], but heuris- 
tically speaking, if the photons have exactly the same mo- 
mentum, they must have infinite uncertainties in their rela- 
tive positions, leading to an infinite pulse width Af . Owing to 
the abrupt increase in 4Af 2 A© 2 when R approaches 1/N, the 
quantum enhancement becomes much more sensitive to loss. 
In the Heisenberg limit of R — > l/N, At — > oo, any loss com- 
pletely detroys the timing accuracy and leads to an infinite 
jitter, according to Eq. 



B. Linear dispersive systems 

If the system is lossy, dispersive, but linear, it is not difficult 
to show that 

Af 2 (z)=Af 2 (0)+C(0) Fdz!P{z') 
Jo 



-A(0 2 (0) 



C(z) 



Aco 2 (z) 



■ C(0) + 2Aco 2 (0) f <feW), 
Jo 

Aco 2 (0). 



(66) 

(67) 
(68) 

The following result can then be obtained from Eq. (l43l after 
some algebra, 

(f 2 (z)) = (f 2 (0)) + (f(0)Cl(0)+Cl(0)f(0)) f/z'Hz) 

f d mz^ 2 > AtHz) 
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(£2 2 (0)) 



N(z 



-(1- 



0- 
(69) 

This result is similar to that in the previous section, except for 
the presence of quantum dispersion and the dispersive spread 
of the pulse width Af (z) that leads to increased jitter. With ini- 
tially coherent-field statistics, (7' 2 (0)) and (n 2 (0)) are given 
by Eqs. (l22l i and (|20| |. respectively, while by similar argu- 
ments, the coherent-field statistics for ^f£2 + £2f) is 

<™ + 6f) coh = ^. (70) 

This leads to the following position variance for a pulse with 
initially coherent-field statistics, 

<^)> coh =^f, (7D 

which still maintains the coherent-field statistics for the dis- 
persed pulse width and the reduced photon number. In the 
high loss limit (az S> 1), the coherent-field statistics is again 
approached regardless of the initial conditions. 

For an initial jointly Gaussian quantum state, on the other 
hand, the normalized version of Eq. d69l is 



R(z) = 



R(0) 



:C 2 



-az 



R(0) 

+ [4Af 2 (0)Aco 2 (0)+4C 2 ] (1 -e- az ), (72) 
where £ is the normalized effective propagation distance, 

t; = Aco 2 (0) r'dz'Piz'), (73) 
Jo 

and4Af 2 (0)Ao 2 (0) is given by Eq. d64b evaluated at z = 0. As 
long as the loss is moderate so that e~ az ^> 1 — e~ az , quantum 
dispersion, given by the term proportional £ 2 /R(0), becomes 
the dominant effect and overwhelms the initial enhancement 
when £ exceeds R(0)/2. 

If the net dispersion f^dz'^iz') is zero, both quantum and 
classical dispersion are eliminated, and the jitter growth be- 
comes identical to that in a non-dispersive and linear system 
given by Eq. doTT ). 
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C. Soliton-like systems 

The previous sections show that coherent-field statistics is 
maintained in a linear system, but as Sec. |IV| clearly shows, 
non-trivial statistics can arise from the quantum dynamics of 
a nonlinear system. The complex evolution of At (z), C(z), and 
Aco(z) in general prevents one from solving Eq. (l43l l explic- 
itly, except for special cases such as solitons. 

If the dispersion is constant and the pulse propagates in the 
fiber as a soliton, C(z) is zero, while At (z) and Aco(z) can be 
regarded as constant if (t 2 ) <C At 2 and (& 2 ) *C Aco 2 through- 
out propagation. Equation (|43l l can then be solved explicitly, 



(f 2 ( z )) = (f 2 (o)> + (a 2 (o)>j3V- 

2j3 2 A« 2 (0) (e az - 1 



Af 2 (0) 
N(0) 



N(0) 



a 1 



z 
a 



(74) 



where (f(O)Q(O) +Q.(0)f (0)) is asssumed to be zero for 
simplicity. If 4(f 2 (0))(i2 2 (0)) = (tt 2 /9)[1/A^ 2 (0)] is also 
assumed for a soliton pulse for simplicity, Eq. d74l i can be 
normalized to give 



R(z)=R(0)e 



81 R(0) 
2% 2 e- az ( e az -l 



:o) \a) 



12 N(0) 



9 A 2 



a- 



z 
a 



(75) 



In the low loss regime with a A <C 1 and az <§C 1, Eq. d75b can 
be further simplified, 



K(z)«*(0) 



81 /?(0) VA 



12 AT(0) 



(OK) 



(76) 



Quantum dispersion is again the dominant effect in this 
regime, while decoherence effects are much smaller, by a fac- 
tor of az approximately. 

Even if the net dispersion is zero and quantum dispersion is 
compensated, the Gordon-Haus effect cannot be completely 
eliminated by dispersion management in the presence of non- 
linearity and may become significant, as the numerical anal- 
ysis in Sec. HVBl shows. An order-of-magnitude estimate of 
Gordon-Haus jitter can be performed by considering soliton 
propagation in a constant negative dispersion fiber, just as 
in the previous case, followed by a dispersion-compensating 
fiber of length L' with positive dispersion coefficient j3'. If 
L' is short, the effective nonlinearity experienced by the pulse 
in the second fiber can be neglected, and A(o{£) can be re- 
garded as constant. Assuming that j3L + j3'L' = 0, the integral 
in Eq. (06]) can be solved to give the Gordon-Haus jitter, 



<f 2 (L + L')) 



GH 



aAco 2 {Q) 
6N(0) 

aA(o 2 {0) 
6N(0) 



j5 2 L 2 (L + L') 
P 2 L\ 



(77) 



The normalized contribution to the squeezing ratio is therefore 

<f 2 (L + Z/)) GH 



(f 2 (L + L')> 



SQL 



(78) 



Compared with the Gordon-Haus jitter at the end of the first 
fiber given by the last term of Eq. 1761 . dispersion manage- 
ment cuts the jitter by half, but the expression maintains its 
functional dependence on the parameters of the first fiber. 
This estimate also justifies the use of Eq. d59l l to estimate the 
Gordon-Haus jitter at the end of the two fibers in Sec. lIV Cl To 
minimize the impact of Gordon-Haus jitter on the quantum- 
enhanced timing accuracy in a dispersion-managed soliton 
system, the condition 



r 3 54 

L <C — 9 



a 



R 



(79) 



is required. 



VI. CONCLUSION 



In conclusion, the decoherence effect by optical loss on 
adiabatic soliton control and on the transmission of quantum- 
enhanced timing information has been extensively studied. It 
is found that an appreciable enhancement can still be achieved 
by the soliton scheme using current technology, despite an in- 
crease of timing jitter due to the presence of loss. It is also 
found that the quantum-enhanced timing accuracy should be 
much lower than the Heisenberg-limited accuracy to avoid in- 
creased sensitivity to photon loss during transmission, and the 
net dispersion in the transmission system should be minimized 
in order to reduce quantum dispersion and the Gordon-Haus 
effect. 

Although the most important pulse propagation effects have 
been considered in this analysis, higher-order effects, such as 
third-order dispersion, self-steepening, and Raman scattering 
lfl4ll might provide further adverse impact on the quantum en- 
hancement if the optical pulse is ultrashort. In particular, the 
inelastic Raman scattering process is expected to be a signifi- 
cant source of decoherence for ultrashort pulses B22I1 . It is be- 
yond the scope of this paper to investigate these higher-order 
effects, but they should be of minor importance for picosec- 
ond pulses and the propagation distances considered in this 
paper. 

Finally, it is worth noting that while this paper focuses on 
optical pulses, the developed formalism is equally valid for 
describing the transverse position and momentum of optical 
beams 12411 and the center-of-mass variables of Bose-Einstein 
condensates [9]. Decoherence by loss of particles in those 
systems can be studied using the formalism developed in this 
paper and parameters specific to those systems. 
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APPENDIX A: DERIVATION OF EXACT HEISENBERG 
LIMITS 

An exact Heisenberg limit can be derived if the inverse 
photon-number operator A^~' is used instead of l/N in the 
definitions of f, Cl, At, and Aco in Eqs. @, (TTB . (TBI , and 
Sl5[ , just as in Refs. and lioll . 



(Al) 
(A2) 

(A3) 
(A4) 



f 




/ rff fA f A, 


&' 






Af' 




y dt t 2 A*A 


Aco' 




J dco co 2 d T 



1/2 



1/2 



These operators are well defined as long as the quantum state 
has zero vacuum-state component ((0|p|0) = 0). Starting 
from the Heisenberg uncertainty principle for T and D.', 



{f' 2 ){Ci' 2 )> 



4 ' 



(A5) 



and the inequality 

' N~ 2 [dco I dco' (CO - Co') 2 a f a' f aa' ) > 0, (A6) 



one can obtain the exact inequality 

(Q! 2 ) <Aco' 2 , (A7) 
and the exact Heisenberg limit for the new position operator, 

(N- 2 ) 



AAco' 2 



(A8) 



The difference between Eqs. ( I28t and (1A81 > is negligible for 
small photon-number fluctuations. The exact Heisenberg limit 
on (iV 2 ) is similar. 



APPENDIX B: NOISE STATISTICS 

In this section the expression (St{z)St(z')) in Eq. d38l is 
calculated. The derivations of (S^(z)Sa(z')) in Eq. d39b and 
(ST{z)Sn(z')+S n (z)S T (z')) in Eq. g0]> are similar. Substi- 
tuting Eq. d35l l into Eq. ( f38l > gives 



(St(z)S t (z')) = j^f jdt jdt' tt' (tfAtfA') + (sUa'V) 



+ <A t M'V) + <A t f S " t A')l, 



where N = N(z), N' = N(z'), § = s(z,t), A = A(z,t), f = 
s(z' ,t'), and A' = A(z',f'). If the noise reservoir is in the vac- 
uum State, s|0 reS ervoir) = (Oreservoirl* = 0, so only the last term 
in Eq. dBU is non-zero, 

(St(z)St(z')) = Jdt Jdt' tt'iPss'U') 
1 

~~ AW 



x 



dt J dt 1 tt' 
(A f A') a8(z-z')8(t - 1') + (AV^A 7 ) 



aAt 1 



8(z-z') 



j^j Jdt J dt' tt' (aV^A') . (B2) 



The first term on the right-hand side of Eq. ( IB2b is the desired 
result, while the second term can be rewritten as 



— - f dt f dt' tt' (AV^A 7 ) 
AW' J J X ; 

= ^JdtJdt'tt'{[A\^) [s,A>]). 



(B3) 



If the system is linear, the commutator between s and A is 
always zero 11311 . but because § does not commute with A* 
and A is coupled to A' by the nonlinear term in Eq. d3"Tl) . $ 
may fail to commute with A. That said, it can be argued that 
the optical field operator must always commute with future 
noise operators due to causality and the infinitesimally short 
memory of s, 



[AV f ] =0ifz<z', 
[§,A'] = if z>z', 



(B4) 
(B5) 



so Eq. ( IB3b can be non-zero only at z — z' ■ The commutator 
between § and A at z — z' due to the parametric coupling of A 
and A* can be estimated by a perturbative technique. Consider 
an integral form of Eq. d3~T| ) with the nonlinear term and the 
Langevin noise term only, 



z+Az 



dz'[iKA\z')A{z')A{z')+s(z')], 

(B6) 



A(z + Az) =A(z) + 
andA'V) 

given by the Hermitian conjugate of Eq. ( 1B6| >, 



(Bl) 



A\z') =A\z) + J dz" [-iKA\z")A\z")A(z")+s\z")\ . 

(B7) 

The commutator between s and A at z + Az becomes 

[f(z + Az),A(z + Az)] 
= iKj' +A ' dz' [S(z + Az) 7 A\z')]A(z')A(z'). (B8) 
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s(z + Az) commutes with A(z') because z + Az > z', while it 
fails to commute with A^(z') because A^(z') given by Eq. jBlj 
depends explicitly on s' . Thus, in the leading order of Az, 

[.f(z + Az),A(z + Az)] 
« iK^ +A ' dz' J dz" [£(z + Az),S^(z!')]A(z')A(z'), (B9) 

which approaches in the limit of Az — > 0. Hence § commutes 
with A at z = z', and the position noise is given only by the first 
term on the right-hand side of Eq. (IB2I ). resulting in Eq. 



APPENDIX C: THE JOINTLY GAUSSIAN STATE 



A Fock state can be expressed as 

I AO = J da>i... J dm <j>{coi,...,a>N)\coi,---,CD N }, 
= J dt x ... j dt N \if(ti,...,t N )\h,...,t N ), (CI) 



where the spectral and temporal eigenstates are given by 

1 



,a>N) 



/Nl 



a T ((Q l )...a T ((Q N )\0), 



(C2) 



|f 1 ,...,f iV ) = ^=A t (f 1 )...A t (f i v)|0). (C3) 

Theses states are eigenstates of the following operators rele- 
vant to our purpose, 

1 



&\ah,...,<a N ) = ( jjY, 03 " )\<*>u---,<*>n), 



(C4) 



f\h,...,t N )= (ij>)|fi,.. .,**>, (C5) 



~ j dco ©^aj ©!,...,©#) 



^E w « )\(Oi,...,CD N ) 



^Jdtt 2 PA\t u ...,t N ) = (if tf)\t u ...,t N )- (C7) 

0(©i, . . . , (On) is the spectral multiphoton probability ampli- 
tude, and it is related to the temporal probability amplitude 
\]/(t\, . . . ,fA?) by the A^-dimensional Fourier transform in the 
slowly-varying envelope regime. Both amplitudes should also 
satisfy normalization and boson symmetry. To study temporal 
quantum enhancement, it is convenient to define the probabil- 
ity amplitude as a jointly Gaussian function [2511 . 

N \ 2 



0(©i, ■ ■ ■ ,(On) = Cexp 
1 

~4P 



j N / j N 

11=1 \ m=l 

N \ 2 



(C8) 



y(h,...,t N ) = Cexp 



^ «— l 



where B and b are arbitrary and real constants, and C and C' 
are normalization constants. Explicit expressions for (i2 2 ), 
(f 2 ), Aco 2 , and Af 2 can be obtained using Eqs. dC4l >-( IC7b and 
Appendix B of Ref. JH, 



(&)=B 2 
(f 2 ) = 



1 



AN 2 B 2 ' 



Aco A 



At 1 



B + \ I 
1 

4N 2 B 2 



A 7 



1 



1 



N Ab 2 ' 



(CIO) 
(Cll) 

(C12) 
(C13) 



In the limit of b — > 0, (r 2 ) reaches the Heisenberg limit, 



<f 2 ) = 



1 



4N 2 A© 2 : 



(C14) 



and the quantum state can be written as a state of photons with 
maximal coincident-frequency correlations, 



\N)<* Jdco exp(-|^r) !©,...,©>. (C15) 



On the other hand, when B 2 = b 2 /N, (f 2 ) is at the standard 
quantum limit, 



(f 2 ) 



1 



ANA® 2 ' 



(C16) 



(C6) the quantum state has only one excited Gaussian mode [25], 



\N) « Jdcoi...J dco N ri ex p( _ ^) •••,<%) 



J dco exp(-^)a f (©) 



|0), 



(C17) 



and therefore also satisfies the coherent-field statistics lfl2l 
1 1311 . These limits and the corresponding quantum states are 
consistent with those suggested in Ref. [1]. With Eqs. (IC12I ) 
and (IC13I ), the pulse width Af can be determined explicitly in 
terms of Aco and the squeezing ratio R = Aco 2 / (NB 2 ), 



At 2 = 



AAco 2 



R (l-l/N) 



2 1 



N l-l/(NR) 



(C18) 
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